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Thermal protection system (TPS) for re-entry vehicle 


Introduction 


qp Quantification of the heat loads on TPS surfaces for re- 
entry vehicles is required to ensure the safety of the crew 

qp TPS surface can exhibit localized roughness, which 
may increase heating due to early transition to turbulence 




Figure: heat shield of the Space Shuttle in orbit during flight STS-114 with gap 
filler sticking out. 
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General configuration 


Physical model and numerical method 


leading-edge shock 



roughness element 
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© Compressible Navier-Stokes equations ! 

© Code: Nagarajan, Lele and Ferziger (2003) [1] 

© High order in space (06, compact) and time (03, explicit) 
o Immersed boundary method to represent roughness 
9 Spanwise periodicity: array of roughness elements 



[1] Nagarajan, S., Lele, S.K. and Ferziger, J.H. (2003), J. Comput. Phys. 191, 392-419. 
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Freestream and forcing parameters 


Physical model and numerical method 


TP Similar setup as in Marxen & laccarino (2008) [1] 
qp Calorically perfect gas, Sutherland’s law with T s 
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TP 3-D roughness (“square”) on an adiabatic wall 
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[ 1 ] Marxen, 0., laccarino, G. (2008), AIAA 2008-4400 
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Steady-state base flow 
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Mean flow I: streamlines & shocks 


Flow over surface with localized 3-D roughness at Ma=4.8 


qp Separation in front < separation in the back 

qp No shock present in x-y planes away from 
the roughness element(s) 


z 



> X 



Figure: Contours ofdp/dx = 0.4 (black) together with contours of u=0 (blue) 
and selected streamlines (grey lines with arrows) in the center plane z = 0. 
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Mean flow II: streamwise vortices 


Flow over surface with localized 3-D roughness at Ma=4.8 


Streamwise vortices (grey) streamwise streak(s) mean-flow gradients (color) 



Figure: isosurface of h 2 (light grey) and recirculation regions (dark grey) 
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Mean flow III: verification 


Flow over surface with localized 3-D roughness at Ma=4.8 


qp Overall good agreement between body-fitted 
(FLUENT) and immersed boundary results 



Figure: Contours of U from body fitted (black contours) and immersed 
boundary method (red contours), planes cutting through the roughness. 
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Mean flow IV: verification (cont’d) 


Flow over surface with localized 3-D roughness at Ma=4.8 


qp Some difference in the spanwise position of the 
streamwise vortices / streaks behind the roughness 
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Figure: Contours of U & T from body fitted (black contours) and immersed 
boundary method (red contours) behind the roughness, x- 18. 
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Mean flow V: transient growth? 


Flow over surface with localized 3-D roughness at Ma=4.8 


qp Streamwise vortices cause a streamwise (u’) streak 


qp Transient growth in individual modes, but no significant 
growth visible in the sum (a non-linear effect?) 
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Figure: streamwise velocity u’ max for F=0 (steady component) for several 
spanwise wave numbers k and their sum. Location of the roughness in grey. 
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DNS of the perturbed flow and comparison with 
stability analysis by Groskopf & Kloker, 2008 [1] 


[1] Groskopf, G., Kloker, M., Marxen, 0. (2008), Proc. of the Summer Program, CTR, Stanford 
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Disturbance evolution I: DNS with 2-D forcing 


Localized 3-d roughness: disturbance flow 


4P 2-D forcing with frequency F=2uf*(jj*/(p*u* 2 ))=0.4U10- 4 
upstream of the roughness (“first mode”) 

4P Fourier analysis in time, disturb, maximum over y&z 



Figure: streamwise velocity u’ max for F=0.41x10' 4 
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Disturbance evolution II: DNS vs. bi-global theory 


Localized 3-d roughness: disturbance flow 


TP DNS vs. bi-global Theory (Groskopf et a!., 2008 [1]) 



Figure: u’ max (max. over z&y) for F=0. 41x1 O' 4 (red: DNS, black: theory) 

[1] Groskopf, G., Kloker, M., Marxen, 0. (2008), Proc. of the Summer Program, CTR, Stanford 
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Disturbance evolution III: amplitude functions 


Flow over surface with localized 3-D roughness at Ma=4.8 


Presence of a y-mode in DNS due to 2-D forcing 
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Conclusions 


Conclusions and outlook 


E 


TP A localized 3-D roughness causes boundary-layer 
separation and (weak) shocks 

TP Most importantly, streamwise vortices occur which 
induce streamwise (low U, high T) streaks 

TP Immersed boundary method (volume force) suitable 
to represent roughness element in DNS 

TP Favorable comparison between bi-global stability 
theory and DNS for a “y-mode” 


TP Outlook: 

- Understand the flow physics (investigate “z-modes” in 
DNS through sinuous spanwise forcing, study origin of the 
beat in DNS) 
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